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ABSTRACT 

We performed K band surface photometry for luminous early-type galaxies in a 
nearby rich cluster ABELL 2199. Combining it with B and R band surface photometry, 
radial variations oi B — R and R — K colours in the galaxies were investigated. It is 
found that the inner regions of the galaxies are redder in both oi B — R and R ~ K 
colours. Comparing the radial variations of both of the colours with predictions of 
Simple Stellar Population (SSP) models for a range of ages and metallicities, it is 
suggested that the cluster ellipticals have negative metallicity gradients but their age 
gradients are consistent with zero, although our sample is small; the typical metallicity 
gradient is estimated to be —0.16 ± 0.09 in dlog Z/dlogr, while the age gradient is 
estimated to be —0.10 ± 0.14 in dlog (age)/(ilogr. Considering that similar results 
have also been derived in the other recent studies using samples of ellipticals in the 
Coma cluster and less dense environments, it seems that there is no strong dependence 
on galaxy environment in radial gradient of stellar population in elliptical galaxy. 

Key v^rords: galaxies: elliptical and lenticular, cD — galaxies: evolution — galax- 
ies: formation — galaxies: clusters: individual: ABELL 2199 — galaxies: fundamental 
parameters (colour gradient). 



1 INTRODUCTION 

Elliptical galaxies in a cluster show a tight colour-magnitude 
(CM) relation; colours of more luminous ellipticals are sys- 
tematically redder (e.g., Bower, Lucey, & Ellis 1992). It has 
been found that the relation still holds in distant rich clus- 
ters around z ~ 1 (e.g., Stanford, Eisenhardt, & Dickinson 
1998). Based on these observational results, galaxy evolu- 
tion models taking into account chemical evolution claim 
that eUiptical galaxies in rich clusters formed at, e.g., z > 3 
through a monolithic collapse of a gas cloud (e.g., Kodama 
& Arimoto 1997). In this monolithic formation, a radial vari- 
ation of stellar population in the sense that stellar metallic- 
ity is higher on average towards the galaxy centre {metal- 
licity gradient hereafter) is expected to form as a result of 
initial starburst and subsequent blow of galactic wind; a 
more extended period of active star formation and thus more 
chemical enrichment is expected in the inner portions of a 
galaxy (e.g., Kawata 1999). Indeed, many elliptical galaxies 
are known to have radial gradients of colours and absorp- 
tion line strengths: colours are redder and metal absorption 
hnes are stronger in the inner regions (e.g., Peletier et al. 
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1990; Davies, Sadler, & Peletier 1993), and the evolution 
of the colour gradients investigated by looking at distant 
ellipticals suggests that these gradients originate from the 
metallicity gradients (Tamura et al. 2000; Tamura & Ohta 
2000; Saglia et al. 2000; La Barbera et al. 2003). 

On the other hand, a bottom-up structure formation 
based on the Cold Dark Matter (CDM) universe gener- 
ally succeeds in reproducing many observational aspects of 
galaxies in the present day universe. In this scheme, an el- 
liptical galaxy is considered to form via a major merger of 
galaxies. Accordingly, even if the progenitors of an ellipti- 
cal galaxy have some radial variations of stellar population, 
the variations are more likely to be smeared out during the 
merger. It is also expected that the histories of mass assem- 
bly and star formation of an elliptical galaxy can be different 
from galaxy to galaxy, even if elliptical galaxies having sim- 
ilar masses at a given redshift are considered. Therefore, the 
radial gradient of stellar population in an elliptical galaxy 
would not be expected to correlate with the galaxy mass, 
although there is a possibility that the gradients in more 
massive ellipticals are less steep because they are likely to 
experience more merging events on average than less mas- 
sive systems. 

These considerations show that investigating colour or 
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absorption line strength gradients in elliptical galaxies in a 
statistical manner is a potentially powerful probe of galaxy 
formation process (Kobayashi 2004). We recently performed 
B and R bands surface photometry for elliptical and SO 
galaxies (E/SOs) in ABELL 2199, which is one of nearby 
rich clusters, to study their colour gradients (Tamura & 
Ohta 2003; hereafter Paper I). From these optical data, it is 
found that (1) the average metallicity gradient in the ellipti- 
cals estimated from their colour gradients on the assumption 
of no age gradient is ~ —0.3 ± 0.1 in dlog Z/dlogr, which 
can be reproduced by a recent model of the monolithic-like 
galaxy formation (e.g., Kawata 1999), and (2) for the galax- 
ies brighter than R = 15 mag [L ~ L* at the distance 
of ABELL 2199), more luminous ones tend to have steeper 
colour gradients (Figure0. These results are consistent with 
the monolithic collapse scenario rather than the bottom- 
up scheme (Larson 1974; Carlberg 1984; Kawata & Gibson 
2003). Considering that almost all the other previous stud- 
ies have targeted ellipticals in environments less dense than 
ABELL 2199 and no such correlations have been found in 
the previous studies, our data may imply an environmental 
dependence of formation process of elliptical galaxy. 

However, it is premature to conclude that the cluster 
E/SOs formed through the monolithic collapse scenario, and 
it needs to be confirmed that the galaxy luminosities cor- 
relate with the metallicity gradients. Although this requires 
to disentangle the age-metallicity degeneracy which is very 
hard in practice, one possible way to address this issue is 
to add another colour information including a near-infrared 
(NIR) band to the optical colour and see the radial varia- 
tions of both of the two colours (e.g., Peletier, Valentijn, & 
Jameson 1990). It will allow us to see whether the colour 
variation is consistent with a pure metallicity gradient or 
a significant contribution of age gradient is essential to ex- 
plain it. In this paper, we present results of K band surface 
photometry for the E/SO galaxies in ABELL 2199 to be com- 
bined with the optical surface photometry (Paper I). K band 
is chosen to make a baseline of wavelength longer to break 
the degeneracy more clearly. Among the luminous ellipti- 
cals plotted in Figure Q we observed 8 galaxies. We briefly 
summarise optical observation and then describe NIR obser- 
vation and data reduction in the next section. Data analysis 
and results are presented in § 3. After discussing the results 
in § 4, we summarise this paper in § 5. 



2 OBSERVATION AND DATA REDUCTION 

2.1 Overall Description of Data Acquisition 

We started studying colour gradients in the cluster ellipticals 
with optical data, which is described in Paper I in detail. We 
performed B and R bands surface photometry for 40 early- 
type (E, E/SO, and SO) galaxies in ABELL 2199 selected 
from the catalog by Lucey et al. (1997). The galaxy mor- 
phology listed by Lucey et al. (1997) is mostly taken from 
Butcher & Oemler (1985) and Rood & Sastry (1972). Al- 
though studying these morphologically selected early-type 
galax;ies could be one option, it is known that there is a 
variety in spheroidal-to-total luminosity ratio among galax- 
ies in each morphological type (Simien & de Vaucouleurs 
1986) and some of them can be disk dominated galaxies. In 



order to securely sample galaxies dominated by spheroidal 
components, we made azimuthally averaged radial surface 
brightness proflles of the galaxies and performed decompo- 
sitions of the profiles into bulge and disk components to 
estimate bulge-to-total luminosity ratios (B/Ts). We sam- 
pled galaxies with B/Ts in B band larger than 0.6. This 
criterion isolates galaxies earlier than E/SO or T ^ —3 ac- 
cording to Simien & de Vaucouleurs (1986), where B/Ts of 
nearby galaxies are also derived in B band. Eventually, 31 
galaxies were sampled out of the 40 galaxies observed. The 
galaxies plotted in Figure are those brighter than R = 15 
mag {L ^ L*) and 8 of these galaxies were observed in K 
band. 

2.2 The optical observation 

Since the details of the optical observation and data are 
described in Paper I, we briefly summarise them here. The 
imaging observations in B and R bands were performed on 

20 and 21 June 2001 with Tektronix 2Kx2K CCD on the 
University of Hawaii 2.2m telescope. One pixel of the CCD 
subtends 0"22 on sky and thus each frame covers a field of 
view of ~ 7' 5 X 7' 5. A typical seeing size during the observing 
run was ~ 0"9 in FWHM. A total exposure time at each field 
is typically 1,800 sec in B band and 750 sec in R band, each 
of which was divided into three or four exposures with the 
telescope dithered. 

The imaging data were reduced with IRAF^ in the stan- 
dard procedure. After bias subtraction and flat-flelding by 
dome-flat frames, sky subtraction was carried out in each 
frame and the sky-subtracted frames at each field were reg- 
istered by sub-pixel shifts. The FWHMs of the point spread 
functions (PSFs) were matched by convolutions of images 
with Gaussian kernels and these images were combined with 
the 3 a clipping algorithm. Finally a PSF size of the stacked 
image in one band was matched to that in another band by 
smoothing with Gaussian kernels. A resulting PSF size is 
1"1 in FWHM. 

2.3 The K band observation 

The K band imaging observation was performed on 20 and 

21 June 2003 with UFTI on UKIRT 3.8m. UFTI has a 
IKxlK array and one pixel of the detector subtends 0"091 
on sky, yielding a field of view of l!5 x 1'5. A typical see- 
ing size during the observing run was ~ 0"7 in FWHM. An 
exposure time of each frame was 40 sec or 60 sec and a to- 
tal on-source integration time is typically 1,200 sec with the 
dither pattern of 9 telescope pointings. Since the apparent 
size of a target is too extended to make fiat-field frames from 
object frames themselves, the on-source exposures were in- 
terlaced by exposures of adjacent blank sky regions with the 
same exposure time as that of the object frames. 

These K band imaging data were reduced with IRAF. 
After dark subtraction, we carried out fiat-fielding using sky- 
fiat frames. To make a sky-flat frame, we stacked ~ 5 blank 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which is operated by the Association of Universities for 
Research in Astronomy, Inc. under cooperative agreement with 
the National Science Foundation 



Optical and near-infrared colour distributions of early-type galaxies in ABELL 2199 3 



J- 

SB 

o -0.05 



-0.1 - 



-0.15 





T 




s ' 1 


71 


J?' 
















6^ 




1 




1 >v 



NGC 6158 
RS 8 
RS 72 
RS 162 
BO 24 
BO 26 
BO 25 
BO 34 



15 14 13 

7? (mag) 

Figure 1. The relationship between B — R colour gradients and R band total magnitudes of luminous {R, < 15 mag or L > L*) early-type 
galaxies in ABELL 2199 (Paper I). Black and grey symbols indicate galaxies with > 3" and those with r'e < 3", respectively. Dashed 
line indicates the trend seen in this plot, which becomes clearer for the galaxies with re ^ 3". The K band data were newly added for 
the numbered galaxies with circles in this plot and these galaxies are discussed in this paper. 
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1.71 
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-0.01 
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-0.03 ±0.01 
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5.9 
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0.04 
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-0.21 ±0.02 


3.0 


0.5 


2.8 
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-0.02 
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-0.04 ± 0.02 


-0.06 ± 0.02 


5.3 


0.3 


1.8 


1.55 


2.78 


0.12 




BO 34 


15.0 


-0.08 ± 0.03 


-0.03 ± 0.09 


6.4 


0.2 


1.3 


1.53 


2.69 


0.12 





Table 1. List of the galaxies with K band photometry. Col. (2), (3), and (5): From Paper I. Effective radii (r^) are estimated 
in the R— band. Col. (6) and (7): The inner and outer cutoff radii are indicated in unit of effective radius. Col. (8) and 
(9): B — R and R — K colours estimated within an aperture of 12" (5 kpc) in diameter, which is equivalent to that in 
Bower et al. (1992) for Coma ellipticals. Col. (10) and (11): The amounts of the offsets applied to B — R and R — K colours, 
respectively. 



sky frames and normalised the stacked image by its modal 
pixel value. An object frame was flat-fielded by the sky-flat 
frame constructed with the blank sky frames which were 
taken close in observation time to the object frame. This 
method allows us to correct for the complex pattern seen in 
the flat-field frame whose contrast slightly varies with time. 
The flat-fielded blank sky frames have very uniform distri- 
bution of brightness across the field of view (global variation 
is ~ 0.1 % or even smaller and pixcl-to-pixcl variation is ~ 1 
%), which is reasonable for an image looking at such a small 
area on the sky. 

Sky subtraction was performed by estimating a single 
background value in each object frame, where the following 
two steps were employed. Firstly, a sky background in each 
frame was estimated by calculating the modal pixel value 
after bright objects were masked and was subtracted. The 
images taken with the telescope dithered around each tar- 
get were registered by sub-pixel shifts and these registered 
images were stacked with the 3 a clipping algorithm. Sec- 
ondly, since objects which are too faint to be found in the 
individual exposure but may cause a wrong estimation of sky 



background are detectable in this stacked image, they were 
picked out using SExtractor version 2.2.2 (Bertin & Arnout 
1996). The coordinates of these objects were de- registered 
to those in the unregistered frames and the pixels where the 
objects should be were masked. The sizes and shapes of the 
masks were defined based on those measured in the stacked 
image. Then, a sky background was re-estimated in each of 
the masked images and it was subtracted. 



These sky-subtracted images were registered by sub- 
pixel shifts and, after the FWHMs of stellar objects in the 
sky-subtracted and registered frames were matched by con- 
volutions of images with Gaussian kernels, they were stacked 
with the 3 a clipping algorithm. Finally, the K band stacked 
image was aligned with the optical data by using stars in the 
field of view and a Gaussian convolution was employed to 
match the FWHMs of stellar objects with those in the op- 
tical images. 
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2.4 Calibration 

For the optical data, the photometric calibration was per- 
formed using the standard stars in Landolt (1992). Sev- 
eral standard star fields were observed at the beginning 
and end of each night. Since the weather was slightly non- 
photometric, an external check of the accuracy of the cal- 
ibration was performed using aperture photometry data 
and growth curves for several galaxies, which were taken 
from HYPERCAT^ There seem to be zero-point offsets be- 
tween our growth curves and those from HYPERCAT; the 
amounts of the offsets are — 0.1 mag in each band. We will 
come back to the treatment of these offsets later. 

For the K band data, the photometric calibration was 
carried out using UKIRT Faint Standard stars (Hawarden 
et al. 2001). About 10 standard stars were observed at the 
beginning and end of each night. The weather was clear and 
stable during the two nights and the photometric zero-point 
can be determined with an accuracy of ~ 0.02 mag from the 
data of standard stars. 

It is found that some of the galaxies have B — R colours 
slightly inconsistent (~ 0.1 mag) with those of typical lu- 
minous ellipticals in the local universe and this can be due 
to the calibration error in B and R bands. We correct for 
these offsets using the CM relation. We assume that the CM 
relation in the Coma cluster also exists in ABELL 2199, and 
their "true" colours are estimated using stellar populations 
modelled by Kodama et al. (1998) using a population syn- 
thesis model so that the CM relations in i7 — 1^ and V — K 
colours obtained in the Coma cluster by Bower et al. (1992) 
are well reproduced. We estimated B — R and R— K colours 
of the sample galaxies within an aperture of 12" (5 /i~^kpc) 
in diameter (Table 0, which is equivalent to that in Bower 
et al. (1992), and estimated the deviations of the colours 
from the CM relations at their luminosities and adopted the 
values as zero-point offsets of the colours. The amount of the 
offset applied to each object is listed in Tabled These are 
corrected for in the following presentations of the colours. 
We note that the R — K colour of BO 24 is ~ 0.2 mag redder 
than the others. Since its B—R colour is not unusual, we also 
attribute this shift to an observational error and corrected 
for it using the CM relation in the same method as above. It 
should be stressed that these corrections are performed by 
shifting zero points of colours and do not affect any spatial 
variations of colour in galaxies, which are focused on in this 
paper. We performed these corrections only for clarity in 
comparing the data with variations of colours predicted by 
population synthesis models as described in the next section. 



3 DATA ANALYSES AND RESULTS 

3.1 Radial Profiles of Surface Brightness and 
Colours in Galaxies 

The radial surface brightness profiles of the galaxies in B, R, 
and K bands are indicated in Figure |5| These profiles were 
obtained along ellipses fitted to the isophotes with the EL- 
LIPSE task in the STSDAS package with a radial sampling 
of 0''22 (1 pix) along major axis. Since bright objects in the 
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neighbor of a target galaxy may have serious effects on the 
fitted ellipses, the neighbors are masked out beforehand. In 
fitting an ellipse to an isophote, ellipticity is set to be a free 
parameter. Galaxy centre is fixed to a centroid of a galaxy 
in the R band image and position angle is also fixed to that 
of an ellipse within which half of the galaxy light is included. 
In the following analyses, all of the radial profiles and the 
parameters related to radius are expressed in terms of the 
equivalent radius of an ellipse: Vah, where a and h are the 
semi-major and semi-minor axes of the ellipse. 

Before obtaining radial colour profiles in the sample 
galaxies, we made the B — R colour maps of the galaxies af- 
ter the seeing sizes in both of B and R bands were matched 
(Paper I). They demonstrate that almost all the galaxies 
show axisymmetric colour distributions independent of mor- 
phology and luminosity and that their radial colour profiles 
and the colour gradients well represent the two dimensional 
colour distributions. We also made the R — K colour maps 
of the galaxies observed in the K band and found that their 
trends are the same as those m B — R. 

The radial profiles of B — R and R — K colour of a 
sample galaxy are constructed by subtracting the surface 
brightness profile in the R band from those in B and K 
bands, both of which were made along the same ellipses fit- 
ted to isophotes in the R band image. The results are shown 
in the bottom panels of Figure |5| An error bar attached to 
each data point in the colour profile includes a local sky sub- 
traction error and a standard deviation of colours along each 
elliptical isophote. To estimate the sky subtraction error, we 
construct a histogram of pixel values in blank regions of a 
fully reduced image where surface photometry of a galaxy 
is performed. Then we adopt the standard deviation around 
the modal value as an error of background estimation. In 
the optical data, the pixel values in an annulus around each 
target whose inner radius and width are ~ 3i?25 and 11" 
(50 pix), respectively, were used to make the histogram. In 
the K band data, pixels which remained after objects were 
masked were investigated in each object frame. In the fol- 
lowing, we focus on the portions of the galaxies between the 
inner and outer cutoff radii, which were defined to avoid 
the regions where colour distributions are seriously affected 
by seeing effects and those where colours are rather poorly 
determined. 

The outer cutoff radius is defined to be a radius 
where the R band surface brightness of a galaxy is 22.5 
mag/arcsec^ (Paper I). At this radius, the typical error of 
B — R colour amounts to ~ 0.1 mag. Since the typical error 
of R — K colour is also ~ 0.1 mag at this radius, it is reason- 
able to adopt the same cutoff as that applied to the optical 
data in investigating the radial profiles oi R — K colour and 
estimating the gradients. The R—K colour profile of BO 
34 (Figurel^J starts becoming noisy at a smaller radius than 
in the other galaxies due to the ~ 40 % shorter integration 
time in the K band. 

The inner cutoff radius is necessary to mitigate seeing 
effects on colour gradients. In order to define this, we per- 
formed simulations using artificial elliptical galaxies added 
on the real images to investigate deviations of measured 
colour gradients from the intrinsic values caused by seeing 
effects. We also investigated their dependences on size, lu- 
minosity, ellipticity, and intrinsic colour gradient of model 
elliptical galaxy (see Paper I for details). Based on the sim- 
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Figure 2. R band image (top), radial profiles of surface brightness (middle), and colour (bottom). In the middle panels, open triangles, 
solid triangles, and circles indicate B, R, and K band surface brightness profiles, respectively. In the bottom panels, triangles and circles 
indicate B — R and R — K colours, respectively, and two vertical lines show the inner and outer cutoff radii. 




Figure 2. Radial profiles of surface brightness and colour - continued. 
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Figure 2. Radial profiles of surface brightness and colour - continued. 



ulations, we defined the inner cutoff radius to be I'.'S, whicli 
can reduce the deviation of a measured colour gradient from 
the intrinsic value down to the level comparable to a fit- 
ting error of the regression line. It should be mentioned that 
since the FWHMs of stellar objects in the K band images 
are smaller than those in the B and R band images and the 
K band data were smoothed to match the FWHMs with 
those in the optical data, this inner cutoff is not required to 
be modified. 

To summarise, the same inner and outer cutoffs defined 
in _B — 7i colour can be adopted to the radial colour profiles 
in R — K. In Table Q the cutoff radii for each galaxy scaled 
by the effective radius are tabulated. The B — R and R~ K 
colour gradients calculated by fitting regression lines to the 
colour profiles between the inner and outer cutoff radii are 
also listed. 



3.2 Spatially Resolved B — R and R ~ K Diagrams 

In Figure |3 the B — R and R — K colours of each target 
galaxy within the region between the inner and outer cut- 
off radii are shown. The colours at outer radii are indicated 
with larger symbols and the radial sampling is the same as 
that in the radial colour profiles in Figure|5|(0!'22). The two 
colour diagrams are ordered in decreasing luminosity from 
left to right and top to bottom (NGC 6158 is the most lu- 
minous and BO 34 is the least luminous). The data points 
in each panel show a basic trend of radial colour variation 
that colours are redder in both oi B — R and R — K at 
smaller radii. It is interesting to note that although the sam- 



ple is small, there seems to be a variety in distribution of 
the colours on the diagram from galaxy to galaxy. 

Although it is hard to disentangle the age-metallicity 
degeneracy, it is still possible to obtain some insight into 
how stellar age and metallicity vary across a galaxy by com- 
paring the B — R and R — K colour variations with predic- 
tions from a simple stellar population (SSP) model. The grid 
superposed on the two colour diagrams indicates variations 
of B — R and R — K colours of an SSP model along with its 
age and metallicity. This is constructed using the SSP mod- 
els by PEGASE Ver 2.0 (Fioc & Rocca-Volmerange 1997) 
with the Salpeter Initial Mass Function (IMF). The details 
of the SSP models for the grid are indicated in Figure|l] We 
sample the ages of 5, 10, and 18 Gyr and the metallicities of 
0.008 (0.4 Zq), 0.02 {Zq), and 0.03 (1.5 Zq). It is stressed 
that only relative values of age and metallicity from a model 
are discussed here and determining their absolute values in 
a galaxy is beyond the scope of this paper. 

Comparing the data points with the grid in Figure O 
it is suggested that the colour distributions tend to be in 
parallel with the prediction for a constant age and it seems 
to be dominated by a metallicity gradient in the sense that 
stars are more metal rich in the inner regions. In most of the 
ellipticals, the data do not seem to support that the primary 
origin of the colour gradients is an age gradient. 

It is worth investigating how much the grid on the two 
colour diagram can vary from model to model. We compare 
the grid by PEGASE with that by Vazdekis et al. (1996; 
hereafter V96) and those by GALAXEV (Bruzual & Char- 
lot 2003) with the Padova 1994 (P94) or Padova 2000 (POO) 
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Figure S. B — R and R — K colour-colour diagrams. The galaxies are arranged in order of decreasing luminosity from left to right and 
top to bottom right (NGC 6158 is the most luminous and BO 34 is the least luminous). The grid superposed with the data points on 
each diagram shows variations of B — R and R — K colours of an SSP along with its age and metallicity predicted using PEGASE for 
the Salpeter IMF. Bold line indicates colour variation predicted for a constant age, and thin line depicts that for a constant metallicity. 
More details of the grid are explained in Figure|l] Note that the same ranges of colours are shown in all panels. 
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Figure 4. Details of the SSP model grid with PEGASE su- 
perposed with the data in Figure |3] The reddening vector for 
Ay = 0.1 mag is also shown for comparison. 



evolutionary tracks in Figure |3 where the Salpeter IMF is 
adopted for all the models. The lower and upper mass cut- 
offs are 0.1 and 100 Mq in PEGASE, 0.0992 Mq and 72 
Mq in V96, and 0.1 Mq and 120 M© in GALAXEV. Note 
that for the GALAXEV model with the P94 evolutionary 
tracks, the metallicity of the most metal rich SSPs is 2.5 
times larger than the solar value, while it is 1.5 times larger 



than the solar value for all the other models. By using PE- 
GASE, variation of the grid due to choice of IMF is also 
investigated. The grids for three representative IMFs are 
indicated in Figure [S] the Salpeter IMF, the Kroupa IMF 
(Kroupa, Tout & Gilmore 1993), and the Millar-Scalo IMF 
(Miller & Scalo 1979). As shown in Figure |3 there is some 
discrepancy in the shape of the grid as well as its location 
on the two colour diagram among the models. On the other 
hand, Figure |S| suggests that both of the shape and location 
of the grid are rather insensitive to choice of IMF. 

Considering contributions of both age and metallicity 
gradients, a colour gradient can be described as follows: 



d{colour) 

dlog r 
where 

_ d{colour) 



^ dlogZ d log (age) 



dloe 



dlogr 



dlogZ 



and 



d{colour) 



(1) 



(2) 



(3) 



a log (age) 

Given the sensitivities oi B — R and R — K colours to age and 
metallicity (i.e., a and h) which can be estimated with an 
SSP model, the age and metallicity gradients can be worked 
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2.5 




^ PEGASE 

GALAXEV(POO)- 



2.5 



1.4 



1.6 

B-R 




Figure 5. Grids constructed with different SSP models are sliown. In each panel, the colour variation predicted with PEGASE (solid 
line) is compared with that predicted with V96 or GALAXEV (dashed line). The Salpeter IMF is adopted in all the models. For the 
grid from GALAXEV with the P94 tracks (right panel), the ages of 5, 10, and 18 Gyr and and the metallicities of Z = 0.008,0.02, and 
0.05 are used. For the other grids, the ranges of age and metallicity are the same as those in Figure |1] 5, 10, and 18 Gyr in age and 
Z = 0.008,0.02, and 0.03 in metallicity. Note that slightly larger ranges of colours than those in Figure l3l are used for the grids to be 
included entirely. 




1.3 1.4 1.5 1.6 1.7 

B-R 




1.3 1.4 1.5 1.6 1.7 

B-R 



Figure 6. Variation of the grid due to choice of IMF is exem- 
plified. In each panel, the colour variation predicted with the 
Salpeter IMF (solid line) is compared with that with the Kroupa 
IMF or the Millar-Scalo IMF (dashed line). All the grids are con- 
structed using PEGASE. 



out simultaneously using the B — R and R — K colour gra- 
dients derived with the data. A similar method has been 
applied to absorption line strength gradients (Kobayashi & 
Arimoto 1999; Mehlert et al. 2003; see also Henry & Worthey 
1999). The coefficient of a is determined within a range of 
metallicity (0.4 < Z/Zq < 1.5) by using colours predicted 
with SSPs having a single age (10 Gyr). The coefficient of 
b is calculated within a range of age between 5 Gyr and 
18 Gyr by using SSPs having the solar metallicity. We cal- 



Model 


B - 


-R 


R- 


- K 




a 


b 


a 


b 


(1) 


(2) 


(3) 


(4) 


(5) 


PEGASE 


0.34 


0.38 


0.82 


0.40 


V96 


0.33 


0.37 


0.63 


0.32 


GALAXEV (P94) 


0.56 


0.39 


1.20 


0.47 


GALAXEV (POO) 


0.27 


0.33 


0.58 


0.41 



Table 2. Sensitivities o{ B — R and R — K colours to age and 
metallicity. The coefficient a (= (d(colour) /d log Z)s,gc) is de- 
termined within a range of metallicity (0.4 < Z/Zq < 2.5 for 
GALAXEV with P94 and 0.4 < Z/Zq < 1.5 for the other 
models) by using SSPs with a 10 Gyr age. The coefficient b 
(= {d{colour)/d log age) z) is determined within a range of age 
between 5 Gyr and 18 Gyr by using SSPs with the solar metal- 
licity. 



culated these coefficients for each of the SSP models (PE- 
GASE, V96, and GALAXEV with P94=' and POO) with the 
Salpeter IMF assumed (Tablel^J and derived age and metal- 
licity gradients of a galaxy. The results from our data are 
shown in Table|3 indicating that negative metallicity gradi- 
ents (stellar metallicity is higher towards the galaxy centre) 
are detected in all the cluster ellipticals but BO 34, while 
an age gradient is significantly detected only in NGC 6158. 
This suggests that the cluster ellipticals have metallicity gra- 
dients while their age gradients are consistent with zero; us- 
ing PEGASE, the typical metallicity gradient is estimated 
to be —0.16 ± 0.09 in dlog Z/dlogr, while the age gradient 
is estimated to be —0.10 ± 0.14 in dlog (age)/dlogr. It is 
mentioned that although estimated values of age and metal- 
licity gradients depend on model, the qualitative trend that 
the cluster ellipticals have negative metallicity gradients but 
age gradients are not detected significantly can be seen in 
all the models. 



^ For GALAXEV with P94, the coefficient a is determined for 
0.4 < Z/Zq < 2.5. 
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Table 3 . Gradients of age and metallicity estimated from the B — R and R — K 
gradients. 



Galaxy 
(1) 


Model 
(2) 


d log (age) /d log r 
(dex-i) 
(3) 


d log Z/dlog r 
(dex-l) 
(4) 


NGC 6158 


PEGASE 
V96 

GALAXEV (P94) 
GALAXEV (POO) 


-0.24 ±0.07 
-0.15 ±0.08 
-0.24 ±0.09 
-0.28 ±0.10 


-0.09 ± 0.05 
-0.19 ±0.06 
-0.05 ±0.04 
-0.09 ±0.09 


RS 8 


PEGASE 
V96 

GALAXEV (P94) 
GALAXEV (POO) 


+0.01 ±0.14 
+0.13 + 0.16 
+0.07 + 0.17 
+0.08 + 0.20 


-0.28 ± 0.09 
-0.42 ± 0.11 
-0.21 ±0.09 
-0.44 ±0.18 


RS 72 


PEGASE 
V96 

GALAXEV (P94) 
GALAXEV (POO) 


-0.12 + 0.14 
-0.03 + 0.16 
-0.10 + 0.17 
-0.11 + 0.20 


-0.16 ±0.09 
-0.27 ±0.11 
-0.11 ± 0.09 
-0.23 ± 0.18 


RS 162 


PEGASE 
V96 

GALAXEV (P94) 
GALAXEV (POO) 


+0.12 + 0.09 
+0.20 + 0.10 
+0.18 + 0.11 
+0.21 + 0.13 


-0.22 ±0.07 
-0.31 ± 0.09 
-0.18 ± 0.06 
-0.37 ±0.13 


BO 24 


PEGASE 
V96 

GALAXEV (P94) 
GALAXEV (POO) 


+0.10 + 0.21 
+0.20 + 0.23 
+0.17 + 0.26 
+0.19 + 0.30 


-0.28 ± 0.14 
-0.41 ± 0.17 
-0.23 ± 0.13 
-0.46 ± 0.26 


BO 26 


PEGASE 
V96 

GALAXEV (P94) 
GALAXEV (POO) 


-0.20 + 0.19 
-0.10 + 0.21 
-0.19 + 0.23 
-0.23 + 0.27 


-0.16 ±0.12 
-0.28 ±0.14 
-0.10 ±0.11 
-0.20 ±0.22 


BO 25 


PEGASE 
V96 

GALAXEV (P94) 
GALAXEV (POO) 


-0.07 + 0.14 
-0.04 + 0.16 
-0.07 + 0.17 
-0.08 + 0.20 


-0.04 ±0.09 
-0.07 ±0.11 
-0.02 ±0.09 
-0.04 ±0.18 


BO 34 


PEGASE 
V96 

GALAXEV (P94) 
GALAXEV (POO) 


-0.33 + 0.33 
-0.33 + 0.39 
-0.39 + 0.43 
-0.46 + 0.49 


+0.13 ± 0.27 
+0.12 ±0.35 
+0.13 ±0.24 
+0.27 + 0.50 



4 DISCUSSION 

It has also been suggested from the other recent studies that 
nearby ellipticals have negative metallicity gradients while 
their age gradients are consistent with zero (Mehlert et al. 
2003; Wu et al. 2004). The presence of a metallicity gradient 
and the absence of an age gradient are consistent with the 
monolithic collapse scenario. The typical metallicity gradi- 
ent estimated from gradients of colours and metal absorption 
line strengths in nearby ellipticals could be reproduced in a 
recent model of galaxy formation through a monolithic-like 
process (e.g., Kawata 1999). Since all stars are presumed to 
form nearly coevally at high redshift in this scenario, no sig- 
nificant variation in stellar age across a galaxy would be ex- 
pected. On the other hand, it is less clear what kind of radial 
variations of stellar populations can be formed in elliptical 
galaxies in the framework of hierarchical galaxy formation, 
while it has recently been suggested by numerical simula- 
tions that metallicity gradients estimated in nearby ellipti- 
cals assuming no age gradients can be reproduced including 



its variety (Kobayashi 2004) . It is intriguing to mention that 
unlike the monolithic collapse scenario, an age gradient can 
be acquired due to star formation associated with merg- 
ers and/or accretions. Exploring what kind of age gradients 
elliptical galaxies could possess in the hierarchical galaxy 
formation may provide clues to understanding their forma- 
tion processes and evolutionary histories (Benson, Ellis & 
Menanteau 2002). 

Our data also suggest that colour distributions on the 
two colour diagram tend to be less scattered in less lumi- 
nous cluster ellipticals. In other words, stellar populations 
may be spatially more uniform in less luminous cluster ellip- 
ticals. This is not because in less luminous ellipticals colours 
are investigated within the smaller portions of galaxies; the 
outer cutoff radii scaled by effective radii are similar in all 
the galaxies (Table 0. If the colour gradients originate from 
pure metallicity gradients, this trend would indicate that the 
metallicity gradients correlate with the galaxy luminosities 
and could support the monolithic collapse scenario (Larson 
1974; Carlberg 1984; Kawata & Gibson 2003). However, the 
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metallicity gradients estimated using both of the B — R and 
R — K colour gradients do not seem to well correlate with 
the luminosities. This may simply be due to the poor statis- 
tics and/or the individual measurements of the metallicity 
gradients may bo too coarse. More precise estimations of age 
and metallicity gradients for a larger sample of galaxies will 
be necessary for detailed studies of the relationship between 
radial variation of stellar population and galaxy luminosity 
in the cluster ellipticals. 

Age and metallicity gradients in nearby ellipticals have 
recently been investigated in a range of galaxy environment. 
In Mehlert et al. (2003), the age and metallicity gradients 
in 35 elliptical and SO galaxies in the Coma cluster were es- 
timated from the gradients of absorption line indices (H/3, 
Mg6 and <Fe>). In Wu et al. (2004), the gradients in 36 
elliptical and SO galaxies in the Sloan Digital Sky Survey 
Early Data Release which are sampled not from a specific 
cluster but from various environments were estimated with 
multi band surface photometry from optical to NIR. Includ- 
ing our study of the luminous early-type galaxies in a rich 
cluster ABELL 2199, the presence of a metallicity gradient 
and the absence of an age gradient have been suggested in all 
of these studies. This may imply that there is no strong de- 
pendence on galaxy environment in radial gradient of stellar 
population in elliptical galaxy. 

Finally, it should be mentioned that dust extinction 
may have some effects on a colour gradient in an ellipti- 
cal galaxy (e.g., Goudfrooij & dc Jong 1995). As a matter of 
fact, even if an elliptical galgixy consists of a mixture of stars 
without any radial gradients of stellar population and dif- 
fusely distributed dust, a calculation of the radiative trans- 
fer in the galaxy suggests that the colour gradients could be 
reproduced only by the dust effects (Witt, Thronson, & Ca- 
puano 1992; Wise & Silva 1996). However, many elliptical 
galaxies show not only colour gradients but also metal ab- 
sorption line strength gradients, which are unlikely to be cre- 
ated by effects of dust extinction. It has also been suggested 
that on average, metallicity gradients in elliptical galaxies 
as estimated by a population synthesis model from colour 
gradients are consistent with those estimated from absorp- 
tion line index gradients (e.g., Pelctier et al. 1990; Davies et 
al. 1993). Nevertheless there are some exceptions and it is 
hard to isolate the effect of dust extinction on an individ- 
ual galaxy basis. Future observation in the far-infrared with 
high spatial resolution may be able to give constraints on 
such a spatial variation of effects of dust extinction. 



5 SUMMARY 

We performed K band surface photometry for luminous 
early-type galaxies in a nearby rich cluster ABELL 2199. 
Combining it with B and R band surface photometry, radial 
variations of B — R and R — K colours in the galaxies were 
investigated. It is found that the inner regions of the galaxies 
are redder in both oi B — R and R — K colours. Comparing 
the radial variations of both of the colours with predictions 
of SSP models for a range of ages and metallicities, it is sug- 
gested that cluster ellipticals have metallicity gradients but 
their age gradients are consistent with zero, although the 
sample is small; the typical metallicity gradient is estimated 
to be — 0.16±0.09 in c?log Z/d\og r, while the age gradient is 



estimated to be — 0.10±0.14 in dlog (agc)/dlogr. Since sim- 
ilar results have also been obtained in the other recent stud- 
ies by investigating ellipticals in the Coma cluster and less 
dense environments, it seems that radial gradients of stellar 
populations in elliptical galaxies and thus their evolutionary 
histories are rather insensitive to galaxy environment. 

Considering the trend found in our optical study that 
less luminous ellipticals have less steep colour gradients, it 
is suggested that they have spatially more uniform distribu- 
tions of both B — R and R — K colours and thus probably 
stellar population. However the metallicity gradients esti- 
mated using both of the B — R and R — K colours do not 
seem to well correlate with the galaxy luminosities. This 
may simply be due to the poor statistics and/or the individ- 
ual measurements of the metallicity gradients may be too 
coarse. More precise estimations of age and metallicity gra/- 
dients for a larger sample of galaxies will be necessary for 
detailed studies of the relationship between radial variation 
of stellar population and galaxy luminosity in the cluster 
ellipticals. 
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